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A method of magnetic resonance perfusion imaging 



FIELD OF THE INVENTION 

The present invention relates to the field of magnetic resonance imaging, and 
more particularly to perfusion tensor imaging. 

5 BACKGROUND AND PRIOR ART 

Diffusion weighted magnetic resonance hnaging is a known prior art 
technique. In dif&ision weighted magnetic resonance imaging the diffusion tensor is obtamed 
&om. the magnetic resonance measurement signal for a defined body region of interest. The 
diffusion tensor is then visualised by known imaging methods (US Patent No. 5,539,3 10; La 

10 Bilhan et al, Diffusion Tensor Imagmg: Concepts and Applications, Jdumal of Medical 
Resonance Imaging 13:534-46 (2001); Lawrence R. Frank, Anlsotropy in High Angular 
Resolution Diffusion- Weighted MRI, Magnetic Resonance m Medicine 45:935-939 (2001),). 
An overview of known diffusion-weighted magnetic resonance techniques is provided by 
Gray L. MacFall J. Overview of diffusion Imaging. Magnetic Resonance Imaging Clm N Am 

15 1998;6:125-138. 

In addition in vivo intravoxel incoherent motion (IVIM) magnetic resonance 
imaging is a known technique (La Bilhan D et al, MR imaging of intravoxel inchoherent 
motions: application to diffusion and perfusion in neurological disorders. Radioloy 
1986:161:401-407.) In particular, the interpretation of FVIM measurements with respect to 

20 classical perfusion has been discussed in the prior art (La Bilhan D et al., Separation of 
diffusion and perfusion in intravoxel incoherent motion MR Imaging. Radiology 
1998:168:497-505; Henkehnan RM. Does IVIM measure classical perfusion? Magn Reson 
Med 1990:16:470-475; Le Bilhan D, Turner R. The capillary network: a link between IVIM 
and classical perfusion. Magn Reson Med Med 1992:27:171-178). This unaging technique 

25 uses gradient encodings for random motion. 

Classical perfusion is a measure of the blood delivered to and used by a 
speciiSed mass of tissue. It is often measured using spin labelling techniques in MRI. In 
contrast, IVIM measures quasi-random blood movement within a single imaging voxel and 
results in a bi-exponential signal attenuation in a standard pulsed gradient spin echo (PGSE) 
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experiment (cf. RJ. Moore et al, in vivo intravoxel incoherent motion measurements in the 
human placenta using echo-planar imaging at 0.5 T, Magnetic Resonance in Medicine, 
43;295-302(2000)). 

It is a common disadvantage of prior measurement techniques that only a 
5 scalar value for the perfusion is obtained even though perfusion usually is anisotropic. A 
scalar perfusion value therefore only gives a limited amount of information on the actual 
nature of the perfusion. 

Therefore the present invention aims to provide an unproved imaging method 
which enables imaging and visualisation of the perfusion anisotropy as well as a 
10 corresponding computer program product and perfusion imaging apparatus. 

SUMMARY OF THE INVENTION 

. The present invention provides for a melliod of perfiision imaging which can 
be used for medical purposes such as abdomuial imaging. 

1 S The present invention is based on the discoveiy that an isotropic perfusion can 

be described by a tensor which is structurally equivalent to a diffusion tensor used for 
dififusion tensor imaging (DTI). In order to provide the data for determination of a perfusion 
tensor a magnetic resonance data acquisition is performed at a low b sensitivity value, 
preferably b=0, as well as at least sbc miagnetic resonance data acquisitions with gradient 

20 encodings in different directions at b sensitivity values of for example below SO [s/mn^], 
preferably between S and IS. 

In accordance with a preferred embodiment of the uivention sbc slopes are 
determined fix)m the 7 measurement points. The slope values form the basis for calculating 
the perfusion tensor, especially its eigenvectors. This is based on the discovery, that for lower 

25 b sensitivity values the signal decay is mostly due to perfusion and not diffusion effects. This 
is why diffusion can be neglected at low b sensitivity values in order to extract the 
. information out of the decay signal which forms the basis for determmmg the shape of the 
perfusion tensor. At least, the sortmg of the eigenvalues and the direction of the eigenvectors 
of the perfusion tensor does not depend much on the diffusion-related signal decay. This 

30 directional information derived from the perfusion tensor can be visualized. For this purpose 
software algorithms which are used in the prior art for diffusion visualization can be reused 
due to the same mathematics governing the visualization of diffusion and perfusion tensors. 

In accordance with a further preferred embodiment of the invention, two 
additional magnetic resonance data acquisitions are performed for higher b sensitivity values. 
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For example one additional magnetic resonance data acquisition is performed for an 

intermediate b sensitivity value of around 200 and a further magnetic resonance data 

acquisition is performed at a high b sensitivity value of around 800. 

For such high b sensitivity values the decay signal is mostly governed by its 
5 diffusion signal component, see e.g. [Petra Miirtz, et al. Abdomen: Diffusion-weighted MR 

Imaging with Pulse-triggered Smgle-Shot Sequences, 258-264, Radiology July 2002]. 

Especially in the abdomen, the diffusion is isotropic, as demonstrated in the cited reference. 

Thus, the diffusion coefficient can be obtained from just two higher b- values measurements. 

For brain imaging it is preferred to perform at least sbc data acquisitions with gradient 
10 encodings for random motion in different directions as diffusion is anisotropic in this case. 

Mathematical analysis also provides the blood fraction value. This fraction is 

known to be isotropic for each measured signal, thus constraining the analysis. Based on the 

diffusion coefficient and the blood fraction value the diffusion component of the decay 

signal can be estimated. By subtracting the diffusion signal component from the 
15 measurement signal Ihe perfusion signal component is obtained. The p^sion coefficient is 

obtained &om the perfusion signal component by determining the signal slope of the 

perfusion signal comjponent, thereto the b=0 and low b-value measurements are analyzed 
In accordance with, a furflier preferred embodunent of the invention only one 

additional magnetic resonance data acquisition is performed for a high b sensitivity value, 
20 leaving out the mtermediary measurement In this case one of the magnetic resonance data 

acquisitions which have been performed for the low b sensitivity values having the highest 

signal value is selected as a replacement for the intermediate measurement. 

It is to be noted that for the mtexmediate and high b data acquisitions no 

gradient encodings in different directions are necessary as diffijsion is isotropic for such high 
25 b values. 

It is a particular advantage of the present invention that due to the limited 
number of data acquisitions, all of the reqmred data acquisitions can be performed for a 
couple of slices through the body (typically 5-10) within a smgle breafli hold, for example in 
less than 16 seconds* 

30 

BRIEF DESCRIPTION OF THE DRAWINGS 

In the following preferred embodiments of the invention will be described m 
greater detail by making reference to the drawings in which: 
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Figure 1 is illustrative of a logarithmic diagram showing magnetic resonance 
imaging signals which form the basis for perfusion tensor determination and unaging. 

Figure 2 is illustmtive of a flow chart for determining the perfusion tensor 
based on the signals shown in figure 1, 
5 Figure 3 is an enlarged view of the decay curves of figure I and is illustrative 

. of a further preferred embodiment of the invention. 

Figure 4 is a blodc diagram of an imaging system. 

DETAILED DESCRIPTION 
10 The general idea of intra-voxel incoherent motion imaging (IVIM) is that the 

observed MR signal decay has a bi-exponential behaviour as function of the diffusion 
weighting factor b, i.e. 

S/So = Al *exp(-b*P)+A2*exp(-b*D), (1) 
where P is the perfusion constant (in mm^/sec, typically 0.05-0.08), and D the diffusion 

15 constant (in mmVsec, typically 0.002). 

It can readily be appreciated that (for h=0) Al+A2=l. Thus, recognizmg that 
the signal related to perfusion can be attributed to (randomly) flowing material in the voxel, 
Al has been designated to represent f, the fraction of flowing material (i.e. blood) in the 
voxel. Hence f denotes the blood fiaction and A2=l-f, which in other words means that the 

20 non-flowing part of the signal contributes to the difKision signal The flection f, being 

representative of the blood content, is an isotropic quantity for each voxel. This knowledge is 
used to constram the data analysis m the rest of this section. 
This leads to 

S/So-f^exp(-bP) + (l-f)*exp(-bD). (2) 
25 The first term of equation 2 will be referred to as the perfusion signal component and the 
second term in equation 2 will be referred to as the difliision signal component. The 
perfusion signal component is decaying much more rapidly than the diffusion signal 
component and thus prfanarily determines the slope of the signal for low b values. 

The perfiision is described by a perfusion tensor, which has the same structure 
30 as the diffusion tensor, i.e. the perfusion tensor is of second rank and symmetric. 

The perfiision tensor can be determined on the basis of the bi-exponential 
signal of equation 2, using the signal decays for low b values whereby the diffusion signal 
component is neglected. This is illustrated by means of figure 1. 
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Figure 1 shows signal decay curves 100, 102 and 104 as a function of b. Magnetic resonance 
image acquisition is performed for b = 0 (point A) as well as for lower b-values (points Bi, 
B2, . . . B6). A good choice for the low b- value is around 10. 

As it can be seen frona the signal decay curves of figure 1 the curves are 

5 approximately linear for low b values. 

For the purpose of determinmg the perfusion tensor the magnetic resonance 
data acquisition for point A is performed without gradient encoding and for the points Bi, B2, 
... B6 at six different gradient directions, respectively. The slopes mi, m2,, . .m^ of the signal 
decay curves for low b-values are determined as illustrated in figure 1 by Imear 

10 approximation. 

By means of the six slope values the perfusion tensor can be calculated by 
means of the same mathematics as used for diffusion tensor calculation (c.f for example 
M.T. Vlaadringerbroek and J.A. den Boen Magnetic Resonance Imaging, Springer- Verlag 
Berlin Heidelberg New York, 1999, in particular section 7.7). 

IS Figure 2 is illustrative of the corresponding flow chart. In step 200 a magnetic 

resonance data acquisition is performed for pomt A (c.f. figure 1) with b = 0 v^thout gradient 
encoding. In step 202 magnetic resonance data acquisition is performed for Bj for a low b 
value with gradient encodmg into a first direction. In step 204 the i is incremented and step 
-202 is performed agam for the next gradient encoding direction. Step 202 is repeated at least 

20 6 times for magnetic resonance data acquisitions into at least six different gradient encodmg 
directions. 

In step 206 &e slope values noj are obtauied fix)m the measurements points A 
and Bi, B2, ... Be by linear approximation, i.e. fi*om the slopes of the straight lines connecting 
A to Bh A to B2 ... A to B6. These slope values provide full mappmg of the perfusion tensor 
25 and are input into the perfusion tensor calculation routine of step 208. 

For increased precision the diffusion signal component of equation (2) is also 
taken mto consideration for the perfusion tensor determmation. One way of accomplishing 
this is illustrated by making reference to figure 3. 

Figure 3 shows the signal decay curves 100 and 102 of figure 1 for higher b- 
30 values. As apparent fi-om figure 3 the signal decay curves become approxunately a straight 
line for higher b values. 

For determination of the diffusion coefficient D two magnetic resonance data 
acquisitions are performed for b values within the linear portion of the signal decay curves. 
For example this can be done for b=200 (point X) and for b^800 (point Y). 
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The slope of the straight Une between points X and Y is the diffusion 
coefficient D. The extrapolation of the straight line between X and Y to b=0 provides 1-f as 
illustrated m figure 3. This way the diffusion signal component of equation (2) is obtamed. It 
is to be noted that the extrapolation by means of the straight line is performed for 
5 convenience of explanation; m a practical application it is preferred to perform the analysis of 
the bi-exponential decay using well Imown mathematical routines such as available firom the 
LAPACK mathematical library. 

The perfusion signal component is analysed by subtracting the estimated 
. diffusion signal component, i.e. Sest=(l-f)*exp(-bD), fi-om the measured value S. For the b=0 
10 value of the curve this provides a value So*f, denoted as So'. For any signal related to non- 
zero b- values, a value S' is derived, A new curve can now be drawn to show ln(SoVS*)- The 
slope of this curve provides P. 

In other words, by subtracting the estimated diffusion signal component, the 
perfusioii signal component is isolated fiom the measured data values; as f is also known the 
IS perfusion coefficient P can be obtained this way. 

For example, this enables a pre-processing of the signal decay curves 100, 
102, 104, ...of figure 1 by subtracting the diffusion signal component before the 
determination of the slopes mi. 

Alternatively the perfusion coefficient P is calculated as explained above for . 
20 every diffiision-sensitisation gradient direction. The resulting set of at least sbc perfusion 
coefiBcients is used to calculate the perfusion tensor map. Specifically, the perfusion tensor's 
eigenvalues, eigenvectors, and rotationally invariant quantities like the trace or the fi^ctiohal 
anisotropy are derived and visualized.. Also, voxels witii substantially co- linear main 
eigenvectors can be connected (starting fix>m a user defined seed position), and the 
25 connecting line can be displayed (tractogn^hy). 

In order to avoid the measurement for point X at the intermediate b value it is 
also possible to use the highest value of the Bi, B2, . . . Be measurements for low b values as 
an approximation. In the example considered here the highest value is Bj (cf. figure 1). In 
order to obtain an approximation for D and 1-^ a straight line between Y and Bi can be used 
3D instead of a straight line between X and Y. 

'One of the advantages of keepmg the number of high b- value acquisitions at a 
minimum is that each acquisition prolongs the breath hold duration by typically 1 to 3 
seconds.Completely avoiding higher b- value acquisitions enables the use of short echo times, 
. because such high b- value acquisitions require much gradient area, thus having a relatively 
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long echo time (typically 60-90 ms, whereas T2 relaxation in tiie abdomen is only 100 ms) . 
The long echo time relative to T2 relaxation significantly lowers the SNR of the 
measurements. Usmg only small b values (<50) significantly reduces the required gradient 
area, and thus enables lower echo times (typically 20-30 ms). This significantly improves the 
5 SNR. 

Figure 4 shows a block diagram of perfiision imaging system 400. Perfusion 
unaging system 400 has magnetic resonance data acquisition device 402. which provides 
magnetic resonance data to computer system 404. For example, the magnetic resonance data 
acquisition is performed by means of a sequence of single-shot echo-planner MR imaging 
10 sequences. 

These sequences are determined by control program 406 which controls MR 
data acquisition device 402 accordingly. The acquired MR data is stored in storage 408. 
Program 410 analyses the acquired MR data which are stored in storage 408 m accordance 
witii the principles as explained above with reference to figures 1 to 3. This way a perfusion 
IS tensor is obtained which is stored in storage 412: The pmiiision tensors stored in storage 412 
are processed by imagmg program 414 which generates a perfusion tensor image which is 
stored in fi^e buffer 416 for display on display unit 418 connected to computer system 404. 
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CLAIMS: 



1 . A method of perfusion imaging comprising: 

- performing a first magnetic resonance data acquisition (A) with gradient encodings for 
random motion at a first sensitivity (b) value, 

- performing a set of at least six second magnetic resonance data acquisitions (Bi, B2, . . . 
S B6) with gradient encodings for random motion in different directions at second 

sensitivity (b) values, 

- determining a perfusion tensor based on the magnetic resonance datai acquisitions. 

2. The method of perfusion imaging of claun 1, the second sensitivity values 
10 being below 50 s/mm^ and the first sensitivity value being substantially snialler than the 

second sensitivity values. 

3. The method of perfusion imaging of claim 1 or 2, whereby the first sensitivity 
value is substantially zero and the second sensitivity values being between five and thirty, 

15 preferably ten. 

4. The method of perfusion hnaging of claims 1, 2 or 3, the magnetic resonance 
data acquisitions being performed by means of a series of single-shot echo-planar magnetic 
resonance data acquisitions. 

20 

5. The method of perfusion imaging of any one of the preceding claims 1 to 4, 
furtiier comprising performing a perfusion tensor visualisation stq>. 

6. The method of perfusion hnaging of claun 5, whereby durectional information 
25 derived &om the perfusion tensor is visualized. 

7. The method of perfusion imaging of any one of the preceding claims 1 to 6, 
further comprising determming of first slope values (mi, m2,. . .nu) between each one of the 
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set of magnetic resonance data acquisitions (Bi, B2, ... Be) and the first magnetic resonance 
data acquisition (A), and determining the perfusion tensor based on the fh-st slope values. 

8. The method of perfusion imaging of any one of the preceding claims 1 to 7, 
5 further comprising: 

- performing of a third magnetic resonance data acquisition (X) at a third sensitivity (b) 
value, 

- performing of a fourth magnetic resonance data acquisition (Y) at a fourth sensitivity (b) 
value, the third sensitivity value being substantially higher than the second sensitivity 

10 values, and the fourth sensitivity value being substantially higher than the third sensitivity 
value, 

- determining of a difiiision coefficient (D) and a fraction value (f) based on the third and 
the fourOi magnetic resonance data acquisitions to provide a difiusion signal component, 

- eliminating of the diffusion signal component fit>m the magnetic resonance data 
IS acquisitions to provide a perfusion signal component, 

- determining of a perfusion tensor from the perfusion signal components. 

9. The method of perfusion unaging of claun 8, whereby a set of at least six third: 
magnetic resonance data acquisitions with gradient encodings for random motion in dififmnt 

20 directions at thh'd sensitivity (b) values is performed, and a set of at least six fourth magnetic 
resonance data acquisitions with gradient encodmgs for random motion in different directions 
at fourth sensitivity (b) values is performed, and the diffusion tensor tensor is determined 
based on the thu-d and the foui^h magnetic resonance data acquisitions to provide a diffusion 
signal component. 

25 

10. The method of perfusion imaging of claim 8 or 9, the third sensitivity value 
being between 100 and 400, and the second sensitivity value being between 600 and 1200. 

1 1 . The method of perfusion imaging of any one of the preceding claims 1 to 7, 
30 further comprising: 

- selectmg one of the second magnetic resonance data acquisitions (Bi, B2, . Be) having 
the strongest measured signal decay, 

- performing a third magnetic resonance data acquisition (Y) at a third sensitivity (b) value, 
the third sensitivity value being substantially higher than the second sensitivity values, 
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- determining of a diffusion coefficient (D) and a fraction value (f) based on the selected 
second and third magnetic resonance data acquisitions to provide a diffusion signal 
component, 

- eliminating the diffusion signal component from the magnetic resonance data acquisitions 
5 to provide a perfusion signal component, 

- determining of a perfusion tensor from the perfusion signal components, 

12. A computer program product, in particular digital storage medium, for 
perfusion imaging comprising program means for determining a perfusion tensor based on a 

10. first magnetic resonance data acquisition (A) and a set of at least sbc second magnetic 
resonance data acquisitions (Bi, 62, ... Be), the first magnetic resonance data acquisition 
being performed at a first sensitivity value and the second magnetic data resonance data 
acquisitions being performed at a second sensitivity value with gradient encodings in 
different directions, whereby tiie first sensitivity value is substantially below the second 

15 sensitivity values, and for performing a perfusion tensor imaging step. 

13. The computer program product of claim 12, the program means being aidapted 
to determine of first slope values (mi, m2,...m6) for each one of the second magnetic 
resonance data acquisitions (Bi, B2, . . . B6) based on the first magnetic resonance data 

20 acquisition (A) and to determme the perfusion tensor based on the first slope values. 

14. The computer program product of clauns 12 or 13, the program means being 
adapted to determine a diffusion coefficient (D) and a fraction value (f) based on thkd (X) 
and fourth (Y) magnetic resonance data acquisitions to provide a diffusion signal component, 

25 to eUminate the difiusion signal component from the magnetic resonance data acquisitions to 
provide a perfusion signal component, and to determme a perfusion tensor from the perfusion 
signal component 

15. The computer program product of claim of claim 14, the program means being 
30 adapted to process a set of at least six third magnetic resonance data acquisitions with 

gradient encodings for random motion in different directions at third sensitivity (b) values, a 
set of at least sbc fourth magnetic resonance data acquisitions with gradient encodings for 
random motion in different directions at fourth sensitivity (b) values, and to determine the 
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difiusion tensor based on the third and the fourth magnetic resonance data acquisitions to 
provide a diffusion signal component 

16. The computer program product of claims 12 or 13, the program means being 

5 adapted to select one of the second magnetic resonance data acquisitions (Bi, Ba, . . . Be) 
having the highest data value, determinmg of a diffusion coefficient (D) and a fraction value 
(f) based on the selected second magnetic resonance data acquisition and a third magnetic 
resonance data acquisition (Y) being perfomied at a third sensitivity value, the thhd 
sensitivity value being substantially above the second sensitivity values, providing a diffusion 
10 signal component based on the diffusion coefficient and the blood fraction value, eluninating 
of the diffusion signal component from the magnetic resonance data acquisitions to provide a 
perfusion signal component, and to determine of a perfusion tensor from the perfusion signal 
components.. 

15 17. A perfusion imaging apparatus comprising: 

- means (402, 404, 406, 408) for performing a first magnetic resonance data acquisition (A) 
at a first sensitivity (b) value and. for performing a set of at least six second magnetic 
resonance data acquisitions (Bi, B2, . . . Be) with gradient encodings m different directions 
at second sensitivity (b) values, 

20 - means (4 1 0, 4 1 2) for determining a perfiision tensor based on the magnetic resonance 
data acquisitions, 

- means (414, 416, 418) for performing perfusion tensor unaging. 
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